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We present a characterization of quantum phase transitions in terms of the the overlap function between two 
ground states obtained for two different vàlues of external parameters. On the examples of the Dicke and XY 
models, we show that the regions of criticality of a system are marked by the extremal points of the overlap 
and functions closely related to it. Further, we discuss the connections between this approach and the Anderson 
orthogonality catastrophe as well as with the dynamical study of the Loschmidt echo for critical systems. 



PACS numbers: 

INTRODUCTION 

Quantum phase transitions (QPT) |1] have drawn a con- 
siderable interest within various fields of physics in the re- 
cent years. They are studied in condensed matter physics be- 
cause they provide valuable information about the novel type 
of finite-temperature states of matter that emerge in the vicin- 
ity of QPT 1 1]. Unlike the ordinary phase transitions, driven 
by thermal fluctuations, QPT occur at zero temperature and 
are driven by purely quantum fluctuations. In the parameter 
space, the points of non-analyticity of the ground state en- 
ergy density are referred to as critical points and define the 
QPT. In these points one typically witnesses the divergence 
of the length associated to the two-point correlation function 
of some relevant quantum field. An alternative way of char- 
acterizing QPT is by the vanishing, in the thermodynamical 
limit, of the energy gap between the ground and the first ex- 
cited state in the critical points. Recently, a huge interest was 
raised in the attempt of characterizating QPT in terms of the 
notions and tools of quantum information |2]. More specifi- 
cally QPT have been studied by analyzing scaling, asymptotic 
behavior and extremal points of various entanglement mea- 
sures Q S 0. 01 • More recently, the connection between 
geomètric Berry phases and QPT in the case of the XY model 
has been also studied |8]. 

In this paper, we aim to provide yet another characteriza- 
tion of the regions of criticality that define QPT. We shall 
show how critical points can be individuated by studying a 
surprisingly simple quantity: the overlap i.e., the scalar prod- 
uct, between two ground states corresponding to two slightly 
different vàlues of the parameters. The physical intuition be- 
hind this approach should be obvious: QPT mark the separa- 
tion between regions of the parameter space which correspond 
to ground states having deeply different structural properties 
e.g., order parameters. This difference is here quantified by 
the simplest Hilbert-space geometrical quantity i.e., the over- 
lap. Note that the square modulus of the overlap is noth- 
ing but the fidelity, widely used in quantum information as 
a function that provides the criterion for distinguishability be- 
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tween quantum states |2]. Therefore, it is a natural candidate 
for a study of macroscopic distinguishability between quan- 
tum states that define different macroscopic states of matter 
(different phases). When applied to cases of many-body sys- 
tems containing many degrees of freedom, the overlap (or, fi- 
delity) might seem to be too coarse quantity, and not bearing 
any apparent information about the difference in order proper- 
ties between quantum phases, to be of any use. Nevertheless 
the main result of this paper is that in some cases it is indeed 
possible to do so. The critical behavior of a system under- 
going QPT is reflected in the geometry of its Hilbert space: 
approaching the QPT the overlap (distance) between neigh- 
boring ground states shows a dramàtic drop (increase). We 
would like to notice that Cejnar et. al. |9| discussed the over- 
lap entropy between the eigen-bases of a system's Hamilto- 
nian and various physically relevant bases, in the context of 
enhanced decoherence effects in the regions of criticality (see 
also El). 

In the following two sections, we conduct our analysis on 
the cases of two simple, yet physically relevant and mathe- 
matically instructive, examples of the Dicke model and the 
XY spin-chain model. Next, we discuss the connection be- 
tween the scaling and asymptotic behaviors and the so-called 
Anderson orthogonality catastrophe lli lli . Moreover, the re- 
lation with the dynamical study of decoherence and quantum 
criticality 1 12] is briefly addressed. Finally, in the last section 
conclusions are discussed. 

For a genèric point in parameter space we use label q E 
K L , where L is the number of external parameters determin- 
ing system's Hamiltonian. As the overlap function depends 
on the difference between parameters as well, we introduce 
q = q + Sq to denote the neighboring point q and the differ- 
ence 5q. Following this notation, we denote the ground states 
by \g) = \g(q)) and \q) = \g(q)). In general, all the func- 
tions F(q) evaluated in the point q we will denote as F, while 
those evaluated in the critical point q c , we will denote as F c 
(note that by combining two cases, we have F c = F(q c + Sq)). 
Then, the overlap function is simply given by the scalar prod- 
uct (g(q) \g{q)) (note that all the results of this paper could be 
easily formulated in terms of fidelity as well). We shall exam- 
ine the behavior of the overlap as a function of q only, while 
keeping Sq fixed and small. 
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THE DICKE MODEL 

Our first example is the Dicke model. It describes a dipole 
interaction between a single bosonic mode à and a collection 
of N two-level àtoms. If for N àtoms we introduce the col- 
lective angular momentum operators J s ,s <E {±,2}, Dicke 
Hamiltonian has the following form (we take H — 1): 



H(X) — ujqJ z + u)à*à + 



X 



•)( 



J+ + J- 



(1) 



Parameter A is the atom-field coupling strength and is the one 
driving the QPT in this model. Therefore, we have q = X 
and denote the Hamiltonian's dependance on that parame- 
ter only. Parameters ujç, and lo stand for the atòmic level- 
splitting and bosonic mode frequency, respectively, while j 
describes the length of a collective spin vector, and is as- 
sumed to be constant and equal to j = N/2. In the ther- 
modynamical limit (N — > 00), which is here equivalent to 
(j — > 00), Dicke Hamiltonian undergoes a quantum phase 
transition for the critical value of its parameter A given by 
A c = (u>ojq)/2. When A < A c , the system is in highly 
unexcited normal phase, while A > À c defines the super- 
radiant phase in which both the field and N àtoms become 
macroscopically excited. The super-radiant phase is charac- 
terized by the broken symmetry given by the parity operator 
II = exp(i7riV), N — (à>à + jz + i): the ground state is dou- 
bly degenerate. As shown in II 311 . by introducing bosonic op- 
erators 6 through Holstein-Primakoff representation fi^ll . the 
above Dicke Hamiltonian Q can be exactly diagonalized in 
the thermodynamical limit. In the normal phase, its form is: 



H 



"(A) = lüoÏÏÍ + ujà^à + X (a f + à) (ÏÏ +b)- juj . (2) 



Its ground state is: g{x,y) = 3 e~ 1/2 < R ' AR > [R = 

(x, y)] where x and y are the real space coordinates as 
sociated to the modes à and b, A = U^MU, M = 

c — s 

diag[e_,e + ] and U an orthogonal matrix U = 

(c = cos 7, s = SÍ117 are given by the squeezing angle 7 = 
(1/2) arctan[4A-y/wa;o / {lo 2 + k>o)])- e ± represent the funda- 
mental collective excitations of the system and are given by: 

4 = \ (lü 2 + u%± ^{lo 2 - uj 2 ) 2 + 16A 2 u; 2 u;^ . From the 

above formula, we see that e_(À c ) = e c _ = 0: the system 
becomes gapless and undergoes a QPT for A = À c . 

The overlap, calculated between two ground states g and g, 
is given by 



{g\g)=2 



\detA detÀ] i 



-=2- 



[deU]i 



[det{A+A)] i [detA] i [det(l+A-U)] * 



(3) 



Note that the overlap is a function of both A and SX. In the 
limit (A — > A c ), with SX > being fixed, det A — s + e- — > 0, 
while det A > det A c = è c + è c L > 0. The same holds 
for det(l + A~ 1 A), for a sufficiently small SX (note that 
lim^A^o A^ 1 = A^ 1 ). But in the present case, it is possible 



to show that for every, and not just small SX, det(l + A^A) 
does not vanish. Using the formula det(l + A) = 1 + 
7ïA + detA for 2 x 2 matrices, we get: det(l + A~ l A) — ► 
1 + Tï(À- 1 A c ) (note that det(À~ 1 A) = [detÀ] -1 detA -> 
0). After a straightforward calculation, we obtain the result: 
Ti^i- 1 A c ) = I || E - [{sc + cs) 2 è% + (ss - ccfè c _] . There- 
fore, Tr(A ( T 1 J 4 c ) > for every A and we can conclude that 
(g\g) (x (e-) 1/4 as (A ->■ A c ). In Ref. |0, it was shown 
that when approaching the critical point from both normal 
and super-radiant sides, the excitation energy e_ drops as the 
square root of A = | A — A c | , which gives us the asymptotic 
behavior of the overlap function in the vicinity of the criti- 
cal point: (g\g) oc A 1 ' 8 . Although we have provided here 
only the results for overlap function for the system in the nor- 
mal phase, the analogous analysis for the super-radiant phase 
gives us the same qualitative results, as the two ground states 
are again the Gaussian-type states, but with translated and re- 
scaled x and y axes. Therefore, we omit it here. 
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FIG. 1: (color online) The overlap function (g\g), equation |3}, as a 
function of A < A c , taken for the resonant case u>o = ^ — 1 and 



5X = 10" 



Note the dramàtic decreasing of the function as we 



approach the point of criticality. 



We conclude this section by presenting the numerical re- 
sults for the overlap function in the normal phase. In Fig. ^ 
we plot the overlap (|3} between two ground states of Dicke 
Hamiltonian for the resonant case ojç, = uj = 1 and SX — 
10~ 6 . We see that it is almost constant and equal to unity for 
wide range of A, apart from the very narrow area around A c , 
when it drastically drops to zero. Such behavior of the over- 
lap function around the point of criticality can be ascribed to 
the fact that the ground state for A = A c becomes completely 
delocalized along one of two rotated axes, as opposed to the 
localized ground state outside of the point of criticality (see 
fil). 



THE XY SPIN CHAIN 

In the following section, we discuss the example of the one- 
dimensional XY anisotropic spin-half chain in the external 
magnètic field. Its Hamiltonian is given by the following ex- 
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pression: 

ff( 7 ,A)=-f; fl±la^^ 1+ ^afay +1 +^) . (4) 

The parameter 7 € R defines the anisotropy, while À € R rep- 
resents external magnètic field along the z axis, up to a factor 
\. Therefore, q — (7, A). The operators trf , a € 
are the usual Pauli operators. This Hamiltonian can be ex- 
actly diagonalized by successively applying Jordan-Wigner, 
Furier and Bogoliubov transformation (see for example (H). 
This way, we obtain the following form of the Hamiltonian: 
H(j,X) = J2kL-M ^k(b{h - 1). The energies of one- 

particle excitations are given by Afe = J e\ + 7 2 sin 2 , 

with e k = cos ^ - A and N = 2M + 1 being the total 
number of sites (spins). One-particle excitations are given by 
the fermionic operators ï>k — cos^dk — is'm ^■(ï i _ k , with 
cos#fc = Efe/A/s. Finally, the ground state |g(7, A)), that 
is defined as the state to be annihilated by each operator bk 
(bk\g(-f, A)) = 0), is given as a tensor product of qubit-like 
states: 

M í 6 9 \ 

|ff(7,A))=®íco8-||0) fc |0)_ fc -ÍBÍn^|l) fc |l)_ fc ) . (5) 

In its space of parameters, the family of Hamiltonians given 
by equation exhibits two regions of criticality, defined by 
the existence of gapless excitations: (i) XX region of critical- 
ity, for 7 = and A G (— 1, 1); (ii) XY region of criticality 
given by the lines A = ± 1 . 

As in the previous example, let us first consider the exact 
overlap function. From equation it follows that the exact 
overlap function between the ground states |<7) and \q) is: 

(g(q)\g(q)) = HcoJ-*^, (6) 

k=i 

where 6k = 0k(q). Note the dependence on the number of 
sites N that is implicit in all the previous formulae from this 
section. In Fig. 0a), we present the numerical result obtained 
using the above equation 0, for N = 10 6 spins and SX — 
5j — 10 -6 . We observe that the regions of criticality are 
clearly marked by a sudden drop of the value of the overlap 
function. As before, we ascribe this type of behavior to a 
dramàtic change in the structure of the ground state of the 
system while undergoing QPT 

In order to investigate the overlap function more quantita- 
tively and relate its behavior to the existence of the regions 
of criticality, we note that while the overlap depends on the 
vàlues of both the parameters q and the difference Sq, the 
regions of criticality are defined by the vàlues of parameters 
only. Therefore, in the following we choose to study the func- 
tions 

k=l V ' k=\ V ' ' 




FIG. 2: (color online) (a) The overlap function (g(q) \g{q)}, as a 
function of A and 7, for TV = 10 6 and 8\ = S-y = 1(T 6 . Note 
the clear dips of the plot in the regions of criticality. (b) S^,(\, 7). 
(c)S£(A, 7 ). 

that define the first non-zero order of the Taylor expansion of 
the overlap function ijóji. Functions 5^r(A,7) and 5^(A, 7) 
are natural candidates for our study because they express the 
"rate of change" of the ground state, taken in the point q. 
They do not depend on the difference Sq, and although for 
every finite N it is possible to find Sq small enough so that the 
exact overlap is arbitrarily well approximated by the expres- 
sion exp(— ^Slj(q)5q 2 ), functions S^(X, 7) and Sj ! (X, 7) on 
their own capture the behavior of the (g(q)\g(q)) function and 
are enough for our current study. They also allow for ana- 
lytic investigation, together with numerical one. In Figs. 0b) 
and|2jc) we present the numerical results for 5^(A,7) and 
Sjf(X, 7), respectively, for N = 10 6 spins. Again, the regions 
of criticality could easily be inferred by simply observing both 
plots. Note that in this case the relative difference between the 
numerical vàlues in the regions of criticality and elsewhere is 
much bigger than in the case of the exact overlap (see Fig. 
|3a)). But now, both plots are needed to detect both regions of 
criticality. This is so because by moving along 7 = 0, while 
keeping |A| < 1, we do not move outside the XX region of 
criticality and therefore do not expect the qualitative change 
in the structure of the ground state, and consequently in the 
behavior of Sjy(A, 7) as well. The same holds for S^(A, 7) 
and the XY region of criticality. 

We first examine the scaling behavior of 5^-(A, 7) and 
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SJf(X, 7) with respect to number of spins N, The numèr- 
ics present us with the following results. First, as expected, 
Sjv(A, 7) and Sjj(X, 7) scale linearly with N when (7 — * 0) 
and (À — > ±1), respectively. In the regions of criticality, we 
have that A| = 1, 7) oc A^/7 2 , while on the other hand, 
for (7 — > 0) we still have Sjf(X, 7) c>c N, without being a 
function of A (we note here that 5^-(|A| < 1, 7 = 0) = 0). 
Such different behavior is a consequence of the fact that while 
the XY region of criticality defines the second order QPT, 
XX is the example of the third order QPT. 

We have also conducted a separate analytical study, con- 
firming the above numerical results. First, we note that for 
every point q in parameter space, and every finite N, partial 

derivatives (7^) and ^7^) are continuous functions of the 
parameters. They can become infinite only in the thermody- 
namical limit, when (N — > 00), and only in the regions of 
criticality. By looking at the explicit form of derivatives (we 
use x k = ^-k): 



7(sinxfc) 



[{cosx k - X) 2 + 7 2 (sina; fe ) 2 



| sin X). | (cos Xk — X) 
[(cosxk - X) 2 + 7 2 (sina; fe ) 2 



(8) 



(9) 



we see that only when the energy (the denominator of 
both of the above expressions) gets arbitrarily small (or zero), 
the derivatives (|8jl and l|9} can become divergent. In other 
words, only when cos Xk gets arbitrarily close to A, and ei- 
ther 7 or s'mxk get close to zero. That is, in the regions of 
criticality. Note that we assume that for every N 6 N and 
k € {1, . . .M}, equation cosxk = A has no solution, which 
presents a genèric case (A's that allow for the solutions of 
this equation form a set of measure zero on the (—1,1) in- 
terval). Therefore, outside the regions of criticality S^(X, 7) 
and Sj í (X, 7) scale linearly with N. 

Regarding the regions of criticality, we first consider the 
scaling behavior of S^(X, 7) in the vicinity of the XX criti- 
cality. As there always exists fco sucn that in the (N — ► 00) 
limit cosxfe — * A, then for such Xk and every finite 7, it 

follows from © that (y^A -> (jsmx^y 1 , when (N -> 

00). In other words, it does not scale with N (note that al- 
though fco = ko(N) is a function of N, limjv^oo sina;^ = 
sin arccos A). As all other derivatives are finite, we have that 
"SjvCM < 1,7 — »• 0) oc N/j 2 . Similar discussion can be 
applied to the case of Sj í (X,-f) in the XY region of criti- 
cality. Again, there exists a qubit defined by k\ = 1 for 
which cos Xk x — * 1 in the (N — > 00) limit, so that its exis- 
tence could bring about the scaling of S^A, 7) larger than 
linear in thermodynamical limit. Using the Taylor expan- 
sion of sine and cosine functions around zero (note that in 
that case, sinxk 1 — > as well), from equation (|9} we ob- 
tain ^Tj^^j oc x^/j 2 — > 0, (N — > 00). In other words, 
1,7) oc N. 



Now, we turn to more interesting cases of the relevant 

functions Sj^(X,j) and S^-(A,7), m the XY and xx re- 
gions of criticality, respectively. Using Taylor expansions 
of sine and cosine functions around zero, we see that in 

' d9 k 



si(\x\ 



X = ±1 the derivative ^-g^H given by k\ = 1 behaves 
like {y^j oc N/(2irj) as (N -> 00) (see equation ©) 

and therefore S^(|A| = 1, 7) oc N 2 /j 2 . We also see that the 
scaling factor depends on 7. Finally, from (|9jl, we also see that 
S^(A,7)<xAT. 

The alternative way to examine the signatures of QPT is to 
look at the asymptotic behavior of two functions near the 
regions of criticality. From the numerical study we obtain that 
the asymptotic behavior of S^(X, 7) in the vicinity of critical 
points A c = ±1, 7 e (0, 1], is given by the following formula: 
S^(A,7) oc 0(7, N)/\l - A| q (t^>. From the study of the 
scaling behavior, we already know that 0,(7, N) — a(j)N 2 
and that 0(7) oc I/7 2 . Further, from numèrics we have that 
the exponent a(j, N) is constant with respect to 7 and ap- 
proaches to a = 1 as (N — > 00). Such asymptotic behavior, 
with constant exponent a = 1 for all 7 £ (0,1] could be 
seen as a consequence of the fact that in that range of param- 
eters the XY model belongs to the same class of universality. 
The numèrics gives also the asymptotic behavior of Sjf(X, 7) 
in the vicinity of 7 = (with |A| < 1) similar to the pre- 
vious one, 5^(A,7) oc b(X, N)/-y^ x ' N \ with the exponent 
/3(A, N) approaching to /3 = 1 as (N — * 00). But, the coeffi- 
cient 6(A, N) depends only on iV, and as noted before, scales 
linearly with it, 6(A, N) oc bN. 



QPT: ORTHOGONALITY CATASTROPHE, LOSCHMIDT 
ECHO 



The above two examples represent a genèric case of a 
many-body system which exhibits continuous QPT only in the 
thermodynamical limit. In the case of the XY model, as the 
number of spins increases, the overlap between two differ- 
ent ground states Q approaches to zero, no matter how small 
the difference in parameters Sq is, so that in thermodynami- 
cal limit each two ground states are mutually orthogonal; they 
live in a continuous tensor product space 1 15]. Such behavior 
of systems having infinitely many degrees of freedom, when 
the two physical states corresponding to two arbitrarily close 
sets of parameters (two arbitrarily similar physical situations) 
become orthogonal to each other, has been already studied in 
many-body physics and is known as Anderson orthogonal- 
ity catastrophe ll 111 . From our study of the XY model, we 
have seen that not only that every two ground states become 
orthogonal in thermodynamical limit, but also the rate of "or- 
thogonalization" between two ground states of large, but finite 
system, changes qualitatively and grows faster in the vicinity 
of the regions of criticality. This way, the regions of criticality 
of an infinite system are already marked by the scaling and 
asymptotic behavior of the relevant functions of a finite-size 
system. Loosely speaking, the regions of criticality of QPT 
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are given as regions where the orthogonality catastrophe is 
expressed on qualitatively greater scale. Notice that recently, 
the occurrence of a particular instance of Anderson-type or- 
thogonality catastrophe was studied for the case of a system 
in the vicinity of QPT 1 16]. 

Now we would like to establish an explicit connection 
between the sort of kinematical approach used in this pa- 
per and the dynamical one of Ref. 11211 . In order to do 
so let us introduce the projected density of states function 
D(u;q,q) = (g(q)\d(u> — H(q))\g(q)) that describes the 
spread of the ground state \g(q)) expressed in the eigenba- 
sis obtained for the point q. Then, the square of the overlap 
can be expressed as 

/>oo 

\(g(q)\g(q))\ 2 = 1 - / D(w)du> (10) 

(Ei denotes the first excited eigenvalue). We see that in re- 
gions in which the spread of \g(q)) with respect to \g(q)) is 
big (quantified in terms of the variance of D(ui)), in other 
words in regions where two ground states differ significantly, 
the overlap will be small. Recently, Quan et al. es- 
tablished a link between the critical behavior of the environ- 
ment and quantum decoherence, showing that the Loschmidt 
echo 1 171 UM L(q, t) of the environment exponentially goes 
to zero as we approach the regions of criticality. A sim- 
ple àlgebra gives us that the Fourier transform of the pro- 
jected density of states is precisely the Loschmidt echo, 
| D(uj)e- luJt duj\ 2 = L{q, t). We see that the kinematics 
of a system, given by the geometry of ground states through 
the overlap function, influences its dynamics as well: the 
smaller the overlap i.e., broader D(ui) the faster the decay of 
the Loschmidt echo. 

CONCLUSIONS AND DISCUSSION 

In this paper, by discussing the examples of the Dicke and 
XY spin-chain models, we have presented a characteriza- 
tion of quantum phase transitions based on the study of the 
scaling and asymptotic behaviors of the overlap between two 
ground states taken in two close points of the parameter space. 
Though this quantity might in general not provide an effi- 
cient numerical tool it is conceptually quite appealing. In fact 
the ground state overlap is a purely Hilbert-space geometri- 
cal quantity, whose investigation does not rely on any a pri- 
ori understanding of the specific kind of order patterns or pe- 
culiar dynamical correlations hidden in the analyzed system. 
While in the case of the Dicke Hamiltonian it was possible to 
analyze the overlap function directly in the thermodynamical 
limit, the case of the XY model is more subtle. The process of 
"orthogonalization" between two ground states has two physi- 



cally different mechanisms in the case when two states belong 
to two different phases: one is a common decrease of the over- 
lap due to infinite number of sub-systems in thermodynamical 
limit, the other is characteristic for the case of QPT and is due 
to different structures of ground states in different quantum 
phases. 
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